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Abstract. - We study resonant transport through a superconducting double barrier structure. At 
each barrier, due to the proximity effect, an incident electron can either reflect as an electron or 
a hole (Andreev reflection). Similarly, transport across the barrier can occur via direct tunneling 
as electrons as well as via the crossed Andreev channel, where a hole is transmitted. In the 
subgap regime, for a symmetric double barrier system (with low transparency for each barrier), 
we find a new T = 1/4 resonance (T is the transmission probability for electrons incident on the 
double barrier structure) due to interference between electron and hole wave- functions between the 
two barriers, in contrast to a normal double barrier system which has the standard transmission 
resonance at T = 1. We also point out as an application that the resonant value of T = 1/4 can 
produce pure spin current through the superconducting double barrier structure. 
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Introduction. — Effects due to the proximity of a 
superconductor have motivated a lot of research in the 
recent past both from theoretical [1-10] as well as experi- 
mental [11, 12] point of view. Due to the proximity effect, 
an electron incident on a normal metal-superconductor 
(NS) interface reflects back as a hole and as a conse- 
c|uence, two electrons arc transferred into the supercon- 
ductor as a Cooper pair. This phenomenon is known as 
Andreev reflection (AR) [13] in the literature of meso- 
scopic superconductivity. An even more intriguing ex- 
ample where the proximity effect manifests itself is the 
phenomenon of crossed Andreev reflection (CAR) [7-9, 14] 
in which an electron from one of the normal leads of a 
normal metal-superconductor-normal metal (NSN) junc- 
tion pairs up with another electron from the other lead 
and as a result, a Cooper pair jumps into the supercon- 
ductor. This nonlocal process can only take place if the 
separation between the points of coupling of the two nor- 
mal metal leads with the superconductor is of the order 
of the size of the Cooper pair itself. From the application 
point of view, the relevance of CAR in the manipulation 
of pure spin currents (SC) [15], spin filter [16] and produc- 
tion of entangled electron pairs in nanodevices [17-19] has 
attactcd a lot of interest in recent times. 

Motivated by these facts, in this letter, we adopt the 
simple-minded definition of SC which is commonly used 



in literature [20-22]. It is just the product of the local 
spin polarization density associated with the electron or 
hole, (a scalar s which is either positive for spin-up or neg- 
ative for spin-down) and its velocity. To generate a pure 
SC in the sense defined above, one can have the two most 
obvious scenarios where (a) there exists an equal and op- 
posite flow of oppositely spin-polarized electrons through 
a channel, such that the net charge current through the 
channel is nullified leaving behind a pure SC , or (6) al- 
ternatively, there exists an equal flow of identically spin 
polarized electrons and holes in the same direction through 
a channel giving rise to pure SC with perfect cancellation 
of the charge current. In this letter, we explore the sec- 
ond possibility for generating resonant pure SC using a 
superconducting double barrier (SDB) structure. 

Proposed device and its theoretical modelling. 

— The configuration we have in mind for the generation 
of resonant pure SC is shown in Fig[T] The idea is to 
design a SDB structure by depositing thin strips of super- 
conducting material on top of a single channel, ballistic, 
one-dimensional (1-D) lead at two places, which can in- 
duce a finite superconducting gap (AiC**^*) in the barrier 
regions as a result of proximity effect of the superconduct- 
ing strips. If the width of the strips is of the order of 
the phase coherence length of the superconductor, then 
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Fig. 1: (a) Cartoon of the SDB structure where a single channel, ballistic, 1-D normal metal (N) lead is connected to a 
ferromagnetic (F) reservoir at its left end and to a normal reservoir at its right. Two patches at the two places on the lead 
depict superconducting material deposited on top of it. (b) Schematic of the potential profile seen by an incident electron. 



both direct electron to electron co-tunneling (CT) as well 
as crossed- Andreev electron to hole tunneling can occur 
across the SDB region. Here it is worth mentioning that we 
restrict ourselves to spin singlet superconductors so that 
both elastic CT and CAR across the junction conserve spin. 
In our theoretical modehng of the system (see Fig. [T|) , we 
first assume that the S-matrix representing the SDB struc- 
ture described above respects parity (left-right symmetry) 
and spin-rotation symmetry, so that we can describe the 
system by an S-matrix with only eight independent pa- 
rameters namely, (i) the normal refection amplitude (tc 
(fc)) for e (h), (ii) the transmission or CT amplitude {tc 
(tc)) for e (h), (Hi) the Andreev reflection (AR) amplitude 
{fAc i^Ac)) for e (h), and (iv) the crossed Andreev reflec- 
tion (CAR) amphtude {tAc (*Ac)) for e (h). (We have used 
the subscript c to denote the composite amplitudes for the 
SDB structure and tilde to distinguish the amplitudes for 
incident holes from the incident electron). If we inject spin 
polarized electron ("f e) from the left lead using a ferromag- 
netic reservoir and tune the system parameters such that 
tc and tAc are equal to each other, it will lead to a pure 
SC flowing to the right lead. This is so because, on an av- 
erage, an equal number of electrons (t e) (direct electron 
to electron CT) and holes (t h) (CAR of electron to hole 
tunnelling) are injected from the left lead to the right lead 
resulting in the cancellation of the average charge current, 
whereas the spins add, giving rise to pure SC in the right 
lead. Note that spin up holes (t h) implies a Fermi sea 
with an absence of spin down electron (which is needed 
for the incident electron (| e) to form a Cooper pair and 
jump into the singlet superconductor). 

Superconducting double barrier. Quantum 
transport in the SDB structure has been studied earlier 
by Morpurgo et al. in Ref. [6]. In their work, they as- 
sumed that both the barriers were reflectionless and also 
that there was no CAR across either of the barriers. They 



then obtained the resultant Andreev reflection and trans- 
mission amplitudes across the SDB by considering multi- 
ple AR processes between the barriers and found a T = 1 
resonance. In this letter, we address the full problem 
allowing all the quantum mechanical processes occuring 
across the two barriers. Hence our set-up is very similar 
to that given in Ref. [6], and comprises of a ballistic nor- 
mal 1-D lead with two short, but finite superconducting 
patches deposited on top of it as shown in Fig[TJ Here 
the structure is connected to ideal ferromagnetic and nor- 
mal electron reservoirs respectively at its two ends. 
and (f>i are the pair potentials and order parameter phases 
on the two patches respectively (i refers to the index of 
the strips). The space dependence of the order parameter 
(which also acts as a scattering potential) for the incident 
electron can be expressed as 

A{x) = A(*)e*e(a;)e(-a; + a) 4- A^'^e*'^^ 

e[x - {a + L)]Q[-x + {2a + L)] (1) 

where, a is the width of the SB , L is the distance between 
the two barriers and is the Heaviside function. 

In contrast to Ref. [6], to obtain the resultant reflec- 
tion, transmission, AR and CAR for the SDB structure, 
one needs to consider all the multiple reflection processes 
in the SDB due to both and vai for each of the barri- 
ers labelled by i. So an electron which enters the region 
between the two barriers has a choice of being reflected as 
an electron or being converted to a hole at each bounce. 
For the numerical analysis, it is more convenient to use the 
alternate method to solve such scattering problems, which 
is to use the standard wave-function matching technique. 

The one-dimensional Bogolubov— de Gennes (BdG) 
equation [23] for (spin up and spin down) electrons and 
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Fig. 2: The variation of Ir^P, l^cP, ^AcP and |tAcP is plotted in units of jh as a function of _E/A for a symmetric SDB system. 
/S./Ef has been chosen to be 0.275. Inset: The region close to one of the resonances is expanded in the inset. Note the existence 
of two closely spaced resonances. 
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The solution of the BdG equations, describing electrons 
and holes with incident energy E inside the normal regions 
(A(,j) = 0), can be written as 
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where, hq^ = yj2m{Ep ± E) and the ± sign in the expo- 
nent of the plane wave solutions corresponds to an excita- 
tion propagating in the ±a: direction. 

Similarly, inside the superconducting barrier regions the 
solutions for electronlike and holelike excitations are 



^±i4 



[x) 



(6) 



(7) 



where, hkf = ^2m{EF ± {E"^ ~ /^.^y^), u± ^^[{l± 

(1 — (A(j)/i?)^)^/^)]^/^, m is the effective mass of the elec- 
tron and Ep is the Fermi energy of the system. 

Hence matching the wavefunctions for the normal and 
superconducting regions (Eq.([3][7])) at the four NS inter- 
faces (a; = 0, a, a -I- i, 2a + L) forming the SDB structure, 
we obtain sixteen linear equations. Numerically solving 
these sixteen equations we obtain the 4x4 §- matrix for 
the SDB structure which, for an incident electron with 
energy E^ can be written as 



tc 

rAc 

tAc 



tc 
tAc 

rAc 



TAc 
tAc 
Tc 
tc 



tAc 

rAc 

tc 



(8) 



In EqlHrc stands for normal reflection of electrons or holes 
and rAc represents AR (reflection of an electron as a hole 
or vice-versa) from the barriers. Similarly, tc represents 
CT or normal transmission amplitude of electrons or holes 
while tAc represents the nonlocal CAR amplitude for elec- 
tron to hole conversion across the SDB structure. The am- 
plitudes depend on the incident energy E, the Fermi en- 
ergy Ep and the length L between the barriers. 

For the normal double barrier system, resonant electron 
transport occurs whenever 6 ~ ±2q'^L = riTr, which is the 
condition for quasi-bound states inside the double barrier. 
The situation for the SDB system is much more subtle. 
Since both electrons and holes are bounced in the normal 
region between the two barriers, there are multiple path 
dependent phases. For instance, an electron that gets re- 
flected as an electron gets a phase of 2q'^L, a hole that 
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gets reflected as a hole gets a phase of 2q~ L, and an elec- 
tron that gets Andreev reflected as a hole gets not only 
the path dependent phase of (q^ — q^)L, but also a A 
dependent phase of cos~^ E/ A [3]. 

We now obtain the resonance condition by using the 
technique of adding up all the Fcynman paths that con- 
tribute to the transmission amplitude. However, here 
since we have both transmission and CAR , we need to 
use matrices for reflection and transmission. Let us as- 
sume the reflection and transmission matrices at each of 
the two superconducting barriers to be the same and given 
by 

R=( ''^ '""^A and T^f/'^ (9) 

\rAeX J \tAe th J 

Note that we have allowed for the amplitudes to be elec- 
tron -electron and hole-hole reflections and transmissions 
re{h) a-nd te(/,) to be different and also similarly electron- 
hole {rAe and tAe) amplitudes to be different from hole- 
electron {vAh and tAh) amplitudes. The A-dependent 
phase X — e^'™** [3] has also been included along 

with each Andreev reflection. The path-dependent phases 
can also be conveniently written in a matrix form as 

where rj = e"^^^ and = e^*'' ^ are the phases picked 
up by the electron and hole respectively, as they move a 
distance L. 

In terms of these matrices, if Ir^l and 0_r,l (each of 
them are column vectors denoting electrons and holes) are 
the incoming and outgoing waves moving towards the left 
or right we find 

(Ol\ ^(R + TPRQPT TQPT \(Ir\ 
\Pr) \ TQPT R + TPRQPTJ {IlJ 

(11) 

where Q = {I - PRPR)'^. We define D ^ Dct(/ - 
(Pi?)^). The condition for resonant transport can now be 
easily found from the composite transmission amplitude. 
For an incident electron going from left to right with wave- 
function given by Ir = (1,0) (II = (0,0)), the amplitude 
for an electron to be transmitted towards the right (upper 
component of Or) is given by the 1-1 component of the 
matrix TQPT = TPQT. The explicit expression for the 
amplitudes arc cumbersome to display, particularly since 
the electron - hole symmetry is broken. Similarly, the 
composite amplitudes for CAR (1-2 component of TQPT), 
reflection (1-1 component of i? + TPRQPT) and AR (1-2 
component oi R + TPRQPT) can also be found from the 
above matrix. 

Clearly, the condition for resonant transport is now 
set by the vanishing of the denominator - i.e.., when 
D{E = Er + iEi) = 0, \tc\'^{Er) has a maximum. Note 
that the composite amplitudes for all the 4 processes have 
the same denominator and hence show resonant behaviour 



when the denominator goes to zero. We find that all four 
of them (using the correct expressions from the matrix) be- 
come 1/4 at resonance, which is a maximum for |tcP, I^AcP 
and |rcp and a minimum for l^^icP- Note also that set- 
ting rA(e,h) ~ tA(e.h) = gi^c thc usual double barrier 
resonance condition, whereas setting r^gjj) = tA(e,h) ~ 
gives the resonance studied by Morpurgo and Beltram [6] . 

Our model docs not explicitly include any exter- 
nal barrier at any of the normal-superconductor or 
superconductor- normal interfaces. In the earlier study of 
SDB without any barrier [6], the approximation A/Ep <C 
1 was also taken, which led to the vanishing of r, and con- 
sequently tA. However, when A/Ep is less than unity, but 
not vanishingly small {e.g., we have taken A/Ep between 
1/4 — 1/10), back-scattering and hence a small non-zero 
value for (which is the reflection at each barrier) does 
exist. We have also checked that the inclusion of normal 
barriers (two external (5-function impurities at each NS in- 
terface) along with thc superconducting barriers does not 
change the result substantially. 

Results. — In this section we describe the conse- 
quences of all the allowed quantum mechanical processes 
across the SDB given by the S-matrix in EqlS] 

Resonance structure:-. As mentioned earlier, we nu- 
merically solve the 16 linear equations obtained by match- 
ing wave-functions at the four NS junctions. We restrict 
ourselves to the subgap regime, where electron energy E is 
much less than the gap energy {E ^ A^ and is small). 
For a symmetric SDB system, Ai = A2 = A; 0i ~ 02. 
The behaviour of |rcp, |tcP, \ tac\'^ and |tAcP as a function 
oi E/ A is shown in Fig. [2l Note that for some particu- 
lar values of E/A, the coherent probabilities for all the 
S-matrix amplitudes given in EqlS] become 1/4. Note also 
that the graphs show two closely spaced resonances. This 
can be understood from the analytic expression in Eqllll 
and more specifically from the denominator D. For small 
values of r, there are two resonances slightly displaced 
from the doubly degenerate pure Andreev level resonances 
which occurs when r = 0. Furthermore, for E Ep, 
(7+ ~ ~ (7^ = \/2mEp, and thc r term in the deter- 
minant for fixed Ep has no significant phase dependence. 
However, the term is multipled by the phase e*''' ~' 
Hence, we have plotted the variation of the transmission 
and CAR probabilities as a function oi 9 = ((7+ — q~)L, 
in FiglSl^a). Note the approximate periodicity of the reso- 
nances is for 9 ^ 9 + 271. 

The width of the resonance depends on the back- 
scattering that occurs at a single barrier, which in turn, 
depends on the value of A/Ep (for r^i) and the scattering 
potential at the barrier (for r,). We have checked numer- 
ically that changing A/Ep between 1/4— 1/10 changes 
rAi which affects the width of the resonances, and chang- 
ing the strength of the (5- function impurities changes r^, 
but otherwise has very little affect on the character of 
the graph. Hence, we have presented our results only for 
A/Ep = 0.275 and A = 0, where A is the strength of the 
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Fig. 3: The behaviour of the spin conductance oc + |iAcP for spin | (j) polarised electrons, is shown, in units of the 
incident spin, as a function oi 6 — (q"*" — q~)L for a symmetric SDB system. The parameter values are given by E/ A = 0.078, 
A/Ef = 0.275 and 01 = 02 = OAtv. The inset shows the variation of \tc\ and |tAc| close to the resonance, in units of e /h. 
Note the closely spaced resonances in the inset. 



(5-function. 

Setting r-i = exactly reverts to the problem studied by 
Morpurgo and Beltram [6] which has only a T = 1 reso- 
nance. This is similar to the usual |tp = 1 resonance of a 
standard normal DB system, (which is obtained by setting 
rAi = 0) although the physical origin is different, since 
for the SDB structure, the electron gets Andreev reflected 
at each boimce between the barriers, instead of getting 
normally reflected. As long as the multiple bounces be- 
tween the barriers occur either through normal reflection 
or Andreev reflection, but not both, the transmission reso- 
nance remains unimodular. On the other hand in our case, 
we have allowed for all the quantum mechanical processes 
that can occur at each barrier. Hence, we have multiple 
bounces between the two barriers involving both r,; and 
rAi- This seems to lead to a new non- unimodular reso- 
nance at T = 1/4, where in fact, all the quantum mechan- 
ical probabilities become 1/4 because all the 4 processes 
show resonant behaviour. This is the main point of this 
letter. It is hence clear that the very occurence of the 
T = 1/4 resonancce requires the presence of all the four 
amplitudes i.e., in a 'single' channel problem with just one 
reflection and one transmission, the only resonance that is 
possible is the standard unimodular resonance. The non 
unimodular T = 1/4 resonance requires the presence of 
two 'channels'. 

Pure spin current:-. As an application of the above 
SDB geometry, we point out that this geometry can be 
used to produce pure SC in a resonant fashion. The pro- 
posal for generating pure spin current using NSN junc- 
tion was discussed earlier in Ref. [15], but there it in- 
volved non resonant production of pure SC unlike the 
present case. In our analysis the spin conductance is de- 



fined as Gs'\{Gsi) oc li^cP + l^cP in units of the inci- 
dent spin, whereas the charge conductance is given by 
Gc^Gci) = {eyh)[\tAc\^ - The T and i arrows 

in the subscript represent the spin polarization of the in- 
jected electrons from the ferromagnetic reservoir as shown 
in Fig[T] The sum of contributions coming from two op- 
positely charged particles (electrons and holes) gives rise 
to the negative sign in the expression for Gci{Gci). The 
interesting point to note here is that for an electron in- 
cident on the barriers, if the amplitudes for the CT and 
CAR are identical, then it will result in equal probability 
for an incident electron to transmit as an electron or as 
a hole across the barriers. This results in the vanishing 
of the charge current. On the other hand, in our geom- 
etry, if the incident electron in the lead is | or ], spin 
polarized, then both the transmitted electron due to tc 
and hole due to tAc will have the same spin polarization. 
This is true because in our analysis we have assumed that 
the superconducting patches are spin singlets and hence 
spin remains conserved. Note, however, that if the in- 
cident electrons were not spin polarised {i.e., if we had 
a normal metal reservoir instead of a ferromagnet), then 
even when CT and CAR are equal, there would be both up 
and down spin electrons and holes transmitted, and hence 
there would be no SC . Therefore, when the symmetric 
SDB structure with a ferromagnetic reservoir is tuned to 

resonance i.e. |tcp = |iAcP = ^ and if a spin polarized 

beam (say "f spin polarized according to FiglT|) of electrons 
is incident on the barriers, then the outcome would be res- 
onant production of outgoing pure SC . In this resonant 
situation 25% of the incident spin-up electrons get trans- 
mitted through the barriers via the CT process and 25% 
get converted to spin-up holes via the CAR process as they 
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pass through the barriers. Hence the transmitted charge 
across the barriers is zero on the average, but there is pure 
SC flowing out of the system. The behaviour of Gs]{Gsi) 
for the SDB system as a function of ((7+ — q')L \s shown 
in FiglJJb). At the resonance, Gs-]{Gsi) becomes 0.5 and 
Gc'\{Gci) becomes for a spin polarized electron beam 
which is a clear manifestation of pure SC in a SDB geom- 
etry. 

Discussions and Conclusions. In this letter, we 
have studied a superconducting double barrier system and 
have shown that one can tune a T = 1/4 resonance in the 
system. It is crucial to have non-zero amplitudes for all 
the four amplitudes, reflection, Andreev reflection, normal 
transmission and crossed Andreev reflection to see this res- 
onance. Note also that we have restricted ourselves to the 
Blonder-Tinkham-Klapwijk approximation [1] of neglect- 
ing the single electron transfer across the barrier. Hence, 
we have shown the resonance only in the thick barrier 
limit, where the transparency is low. A similar resonance 
was already noted by some of us [24] in a stub geome- 
try where all these four amplitudes were non-zero. The 
special value of T = 1/4 was also noted earlier by some 
of us [10, 15] in the context of a weak interaction renor- 
malisation group study of a NSN junction, where a non- 
trivial fixed point was found which had T 1/4. In all 
these contexts, not only does T have a non-unimodular 
value, all the other amplitudes also have a value of 1/4, 
(jiAcj^ ~ |?'c|^ ~ Vac]"^ = 1/4) as required by unitarity. In 
all these contexts, pure SC is the outcome (at resonance or 
at the fixed point), since the charge current gets nullified 
on the average. 

As far as the practical realization of such a SDB struc- 
ture is concerned, it should be possible to fabricate such 
a geometry by depositing thin strips of a spin singlet su- 
perconductor (like Nb with A ^^l.SmeV [11]) on top of a 
ballistic quantum wire (with Ep ~leV [25]) or a carbon 
nanotube at two places. The width of the strips should 
be of the order of the superconducting phase coherence 
length (10— 15nm in case Nb). The T = 1/4 resonance 
in this SDB geometry can be tuned by varying the energy 
of the incident electron (which can be done by applying 
a small bias voltage between the two reservoirs keeping 
within linear response, so that our calculations are valid) 
for fixed distance {'^ O.bfim) between the two barriers or 
the distance between the two barriers for fixed incident 
energy. However, inclusion of electron— electron interac- 
tion, finite temperature, finite bias, etc can lead to very 
interesting physics in the presence of resonances which is 
beyond the scope of the present work. 

In conclusion, in this letter we have studied resonant 
transport through a SDB stucture where, at an energy 
scale much below the superconducting gap A, probablitics 
for all the coherent amplitudes become 1 /4 in the tunnel- 
ing (or thick barrier) approximation. As an application 
we have also discussed the possibility for production of 
resonant pure SC in this structure. 
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